Charge screening and magnetic anisotropy in metallic rare-earth systems 
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The calculation of magnetic anisotropy constants is per- 
formed beyond the point charge model for a continuous charge 
density distribution of screening conduction electrons. An im- 
portant role of the non-uniform electron density, in particular, 
of the Friedel oscillations, in the formation of crystal field is 
demonstrated. Such effects can modify strongly the effective 
ion (impurity) charge and even change its sign. This enables 
one to justify the anion model, which is often used at dis- 
cussing experimental data on hydrogen-containing systems. 
Possible applications to the pure rare-earth metals and RC05 
compounds are discussed. The deformation of magnetic struc- 
ture near the interstitial positive muon owing to the strong 
local anisotropy, and the corresponding contribution to the 
dipole field at the muon are considered. 



The old problem of strong magnetic anisotropy in the 
rare-earth-based intermetallic systems [jlj, which has a 
great practical importance, is up to now extensively in- 
vestigated. Recently the anisotropy induced by intersti- 
tial hydrogen atoms ^,|| and positive muons Qj has been 
discussed. 

It is accepted now that the main mechanism of mag- 
netic anisotropy origin in the rare-earth systems is the 
crystal field one. Estimations of the anisotropy constants 
are usually performed in the point-charge model. At the 
same time, this model leads frequently to difficulties and 
contradictions with experimental data. For example, the 
calculated anisotropy constant K\ in RC05 compounds 
turns out to be very large and have an incorrect sign [^) . 
Thus screening of the crystal field should play an impor- 
tant role. This screening is often taken into account by 
introducing the effective ion charge Q* which can differ 
considerably from the bare ion charge. In particular, the 
fitted charge of hydrogen ions for the RCo2-based systems 
turns out to be negative (the anion model (HJD). Physi- 
cal explanation of such situations is not simple. Shield- 
ing of crystal fields in ionic solids was treated in Ref. M . 
Screening by conduction electrons, which form some lo- 
calized levels near the rare earth ion, was discussed in 
the paper g. However, these effects turn out to be in- 
sufficient for explaining experimental data. 

In the present paper we perform calculations of mag- 
netic anisotropy constants with account of continuous 
charge-density distribution in a metal. We shall demon- 
strate that the magnetic anisotropy is strongly influenced 
not only by the total charge induced by surrounding ions 
or impurity centres, but also by a concrete form of screen- 
ing electron density, in particular, by the Friedel oscilla- 
tions. 



We consider the magnetic ion at the point r = in the 
crystal field of the surrounding charges. The spherically 
symmetric potential V c f, which is induced by a center at 
the point r = R and originates from the point charge Qq 
and conduction electron charge density Z(\r — R|), has 
the form 



Vcf(T) 



Qo + Qei(\r-H\) 
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where 



Q el (r) = 4vr f 
Jq 



p 2 dpZ(p) 



(1) 



(2) 



is the conduction electron charge inside the sphere with 
the radius r, Q' e i(r) = A-Kr 2 Z(r), the system of units 
where the electron charge e = — 1 is used. At large r 
complete screening take place, so that Q e /(oo) = —Qo- 

The anisotropy constants are determined from the an- 
gle dependence of the energy of the magnetic ions in the 
crystal field 
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cf 



-K x cos 2 o 



(3) 



First we discuss the point charge model. The dependence 
(^) is calculated with the use of the expansion in r, which 
corresponds to the expansion in small radius of /-shell: 
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2K Z 
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Here is the angle between the vectors R and r, which 
can be expressed in the spherical coordinate system of 
the crystal as 



cos = cos cos Or — sm sin Or cos( 



to.) (5) 



where Or and <j)R are the polar and azimuthal angles of 
the vector R. Substituting (||) into (0) and picking out 
the term, which is proportional to cos rO, one obtains for 
the point-charge contribution (cf. Refs. (?]||) 



Kf = -3A(r 2 )ajJ(J-l/2)Q e 2 



(6) 



where (r 2 ) is the average square of the /-shell radius, J 
is the total angular momentum of rare-earth ions, aj is 
the Stevens factor, 



R 



3 cos 2 R - 1 
R~ 3 



(7) 
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For the hep lattice with the parameters c and a (pure rare 
earth metals) we derive after summation over nearest- 
neighbor magnetic ions the standard result (see Ref . || ) 



A = 6a" 



1G 



y 2 ~ 2/3 

(4/3 + y 2 )5/ 2 



2.44a~ 3 ( v /873- 2 /) 
(8) 



where y — c/a ~ 1.57 -j- 1.59 for heavy rare earths Q. 
The small geometrical factor y/8/3—y = 1.633 — y ~ 0.05 
occurs due to that the contributions from the magnetic 
ions in the same plane and in neighbor planes have op- 
posite signs and almost cancel each other. 

For the RC05 compounds, where R = a/v3 for six 
neighbor Co ions in the same plane and R — \\j 'a 2 + c 2 
for six Co ions in the upper and lower planes, we have j^] 



A = 6c 
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2y 2 



(l+?/2)5/2 



3 3/2 



23.4 



(9) 



with y — c/a ~ 0.8, a ~ 5A. In this case the cancellation 
of contributions from different planes does not play a 
crucial role, so that the expression in the brackets is never 
small. 

Consider the case where anisotropy is induced by ran- 
dom next-neighbor substitutional or interstitial impuri- 
ties. For an impurity in the octahedral interstitial of the 
hep lattice we have R = a/y/2, cos#r = l/y/3. Then 
A = for the ideal hep lattice, and the impurity contribu- 
tion to Ki vanishes, as well as (||) (note that the situation 
is different in the hydrogen-containing RC02 systems M ) . 
However, there occurs the local anisotropy term 



S£ c f = -ifllocCOS 2 



(10) 
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Then we derive 

Ki/K? = Q*/Q (14) 
with the effective charge 

Q* = Qo + Qei(R) - ^ttR 3 [Z(R) - RZ'(R)} (15) 

Thus we have obtained the expression for the effective 
charge Q*, which should be used at calculating the ob- 
servable anisotropy constant K\. We see that the ratio 
(|14|) depends explicitly, besides the total charge inside the 
sphere with the radius R, also on the charge density Z(R) 
and the derivative Z'(R). The latter quantity can be large 
provided that the rare-earth ion lies in the region where 
the charge density changes sharply. Note that for the con- 
stant charge density Z (Z'(R) = 0, Q e i(R) = t-ttR 3 Z) we 
have Ki = Kf c and screening is absent (this is connected 
with that the charge in the sphere increases as R 3 , but 
the quadrupole interaction decreases as i?~ 3 ). 

To obtain the value of Q* , one has to investigate 
the charge screening for a concrete electronic spectrum, 
which is, generally speaking, a very difficult task. We 
discuss the one-centre screening problem within a simple 
model of free conduction electrons [E = k 2 /2) in the 
impurity-induced rectangular potential well which has 
the width d and depth E — k 2 /2 [|| . This model enables 
one to calculate the charge distribution of screening con- 
duction electrons in terms of the scattering phase shifts 
r/;. The value of ko should be determined for given kp 
and d from the Friedel sum rule 



K p n c oc = l2V2(r 2 f )ajJ(J - l/2)Q e 2 /a 3 . (11) 

Note that for the positive charges Qo the signs of -fTnoc 
and K\ coincide. 

Now we treat the case of a continuous charge distri- 
bution of the screening conduction electrons. Similar to 
(|), the potential ([[]) can be also expanded at small r. 
Performing the expansion of the integral (||) in 

|r-R| -R = -rcos8+ (r 2 /2R) sin 2 9 + ... 
we obtain up to r 2 

Q el {\r - R|) = Q el (R) - 4nR 2 Z(R)r cos 6 + 2txRt 2 

x {Z{R) + [Z{R) + RZ'{R)] cos 2 6} (12) 

Taking into account the expansion (|J) we have 

Q d (|r-R|) Q el (R) 



Rl 



-R- [1 -R COa0+ 2R^ {3C ° s2 ' 



1)] 



-2Tr[Z(R)(2rRcos 6 - r 2 sin 2 0) - RZ\R)r 2 cos 2 9} (13) 



(16) 



The phase shifts are calculated as 



kdn\(kd) — Pono(kd) 

tan?/ ;>0 (fc) = fcrf J'- 1 ( fcrf ) I Mi(kd) ^ 
kdni-\(kd) — /3ini(kd) ' 



kd 



jo(kd) 
Ji-i(kd) 
ji(kd) 



(18) 



where ji{x) and ni(x) are the spherical Bessel and Neu- 
mann functions, k 2 — k 2 + k 2 . The disturbance of the 
charge density is given by 

1* 00 

5Z{r) = - kdkSp{k,r) 
Jo 

For r > d one obtains 
Sp(k,r)/p (k) = 5>Z + \){[n 2 {kr) - j 2 (kr)]sin 2 m 



(19) 



-ji(kr)ni(kr) sin2r/ ; } 



(20) 



with po(fc) = k/ir 2 . For r < d, joining of the wavefunc- 
tions at the boundary of the potential well r = d yields 
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Sp(k, r)/po (k)= ( 2l + Wi (kd) cos m 
i 

-ni(kd)sm m ] 2 \ji(kr)/ji(kd)] 2 - 1 (21) 

The parameter d should be determined by the geom- 
etry of the lattice near the impurity. In Ref . || , where 
impurities in the Ag host were considered, d was chosen 
to be equal to the Wigner-Seitz radius, so that k F d = 2. 
In the case where the impurity is localized in an intersti- 
tial the choice of d may be different. 

We have performed the calculations for kp d — 2 and 
k F d = 3. At Qq = 1 Eq.(|l|) yields k d = 1.46 and k d = 
1.235 respectively. With increasing kp, the number of I 
values to be taken into account in the sums grows; for 
k F d = 3 the contributions up to I = 3 are appreciable. 
The results are presented in Figs. 1-3. Note that a weak 
singularity occurs at the joining point r = d. 

One can see that at R < d, except for the case of very 
small R where Q*(R) slightly decreases, the derivative 
term results in that Q*(R) grows (despite an increase 
of \Q e i(R)\)- For R ~ d, where Z' is maximum, the 
non-uniform distribution of electron density leads to that 
the effective ion charge Q* is positive and exceeds con- 
siderably its bare value Qo (Qo = 1 hi our case). At 
the same time, with further increasing R the situation 
changes drastically: Z' decreases and becomes negative, 
so that "overscreening" of the ion charge occurs. Such a 
situation corresponds to the "anion model" . In the sim- 
ple model under consideration, the values about —0.5 can 
be obtained for Q* (the value of -1 was assumed for the 
hydrogen ions in Ref. ||). At large distances Q* tends to 
zero, but considerable oscillations of the effective charge 
sign take place, which attenuate rather slowly. One can 
see from Fig. 2 that the contribution of the derivative 
term in ( |i~5| ) predominates at large R too. Of course, 
more complicated models are required for a regular lat- 
tice of screened charges, but the conclusion about an im- 
portant role of the non-uniform distribution of electron 
charge should hold also in that case. 

The effects of screening under consideration can be 
important for the values of the anisotropy constants in 
both stoichiometric rare-earth compounds and systems 
containing hydrogen and other impurities. The experi- 
mental data on the first magnetic anisotropy constant of 
the heavy rare-earth metals and corresponding RC05+Z 
compounds at low temperatures are presented in Table 
1. The values of the effective charge Q* are calculated by 
using (fl4"|), (||). One can see that the values of Q* for the 
pure rare earths are of order of unity. At the same time, 
the effective charges of the Co ions in the RCos+a, systems 
turn out to be negative and very small (a similar situa- 
tion takes place in the case where R is a light rare earth 
Q), so that a weak "overscreening" occurs. Note that in 
fact the effective charges may be different for two kinds of 
the Co ions which contribute (^), since the corresponding 
values of R are different. Besides that, a possible role of 
the neighbor rare-earth ions in the crystal field formation 
should be considered in such a situation. 



Table 1. The total angular momenta J, Stevens fac- 
tors aj, average squares of the /-shell radius (r 2 ) (atomic 
units) and c/a ratios for rare-earth elements; the experi- 
mental values of K\ (K/rare-earth ion) for the pure heavy 
rare earths [|| and R Cos +x systems (the contribution of 
rare-earth sublattice |l|]), and the corresponding calcu- 
lated values of Q* . The values of K^ xp for Er and Tm are 
obtained from the anisotropy of magnetic susceptibility 
(see Ref. ]l2|] ). The values of the local anisotropy con- 
stant -Kiioc for an impurity in the octahedral interstitial 
of a rare-earth metal are calculated for Q* — 1. 
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Another interestin g e xample is the case of imbed- 
ded positive muons |l^,|l3| which induce a strong lo- 
cal anisotropy. This leads to strong deformation of the 
magnetic structure near the muon and to a considerable 
contribution to the dipole field at the point of its loca- 
tion QQ. According to Table 1, for Q* ~ 1 the local 
anisotropy constants make up about 10 3 K and are large 
in comparison with the values of the molecular field act- 
ing on the rare-earth ions, A ~ 100K. Therefore a strong 
cast of magnetic moments should take place in the limit 
of large Ku oc . In the case K\ > 0, A'n oc > for the 
ferromagnetic ordering along the z-axis we have for the 
dipole field 

AB| ip = AVSM/a 3 ~ 25 kG (22) 

where M ~ y/ J (J + 1)/i_b is the magnetic moment of 
rare-earth ion. In the case K\ < 0, Ki\ oc < the direc- 
tion and value of the dipole field depend on the relation 
between K% and A (see Ref. Q), but lAB^ipl has the 
same order of magnitude. Thus the local distortion of 
the magnetic structure yields the contribution to ABdi P , 
which is of the same order as the sum over the regular 
lattice (cf. the calculation for holmium |^5); note that 
the nearest-neighbor ions give zero contribution to this 
sum in the absence of the distortion). Phenomena which 
occur in ferromagnetic and helical phases at intermedi- 
ate values of Ku oc with changing the relation between 
if Hoc, Ki and A (e.g., with the change of temperature or 
external magnetic field) are discussed in Ref. Q. 

To conclude, the effects of screening of ion (impu- 
rity) charge by conduction electrons turn out to be very 
strong: the anisotropy constants can be strongly modi- 
fied and even change their signs. More quantitative in- 
vestigations with account of a real electronic structure, 
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multi-center effects and non-spherical charge-density dis- 
tribution seem to be important for the problem of mag- 
netic anisotropy in the systems under consideration. 

The research described was supported in part by the 
Grant N 96-02- 16999- a from the Russian Basic Research 
Foundation. 

Figure captions 

Fig.l. The distance dependence of the function 
4ttR 3 5Z(R) for k F d = 2 (solid line) and k F d = 3 (dashed 
line) . 

Fig. 2. The dependence 4:TtR 4 SZ' (R) for the same pa- 
rameter values as in Fig.l. 

Fig. 3. The distance dependence of the effective charge 
Q*(R) for the same parameter values as in Fig.l. 
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